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New York Blower 
Graphs Fan Performance 

by Terry Davis 

TEKniques staff 

with Tim Barnhart 

The New York Blower Co. 

LaPorte, IND 



Designing and building a system to move 
air is no simple task, whether the system 
is for air conditioning, fume evacuation, 
combustion air, or other uses. There are 
many factors to take into account in the 
design, several that relate directly to the 
fan to be used in the system: torque, 
horsepower required, and air flow vol- 
ume and pressure are some of these. 
How, for instance, does the horsepower 
requirement relate to the delivered flow? 
Or how does air flow relate to static 
pressure? And how do these relation- 
ships change over the operating range of 
the fan? Fitting the correct fan into the 
system requires integrating these factors 
into the system design. 

But that leads to another question: how 
does the fan maker communicate these 
parameters and their changing relation- 
ships to the air system designer? At The 
New York Blower Company, graphics 
were chosen as the answer some years 
ago. Fan Performance Curves and 
Speed-Torque Charts were developed to 
be presented to inquiring customers, de- 
pending upon their specific needs. At 
first these charts were calculated and 
plotted by hand. But starting four years 
ago, the 4050 Series Graphic Computing 
Systems began to take on that burden, 
with the help of the 4662 Plotter. 

The New York Blower Company is located 
in LaPorte, Indiana, near the shores of 
Lake Michigan. From this quiet midwest 
setting, they manufacture and distribute 
industrial blowers throughout the United 
States and, through licensees, interna- 
tionally as well. Blower sizes range from 
the smallest, an eight-inch wheel, to their 
huge 85" wheel blower. Three basic cen- 
trifugal fan types are produced — for- 
ward curve (or squirrel cage), radial 
blade wheels, and curved radial blade 
wheels. Other special-purpose fans pro- 



duced include vaneaxial fans, pressure 
blowers, tubular acoustafoil fans, and 
corrosion-resistant fiberglass fans for 
fume exhaust applications. 



Basic Graphic Applications 

There are two 4050 Series Graphic Com- 
puting Systems at New York Blower 
Company. The first, a 4051, has been 
there for about four years, and currently 
resides in the Research Laboratory. The 
recently-acquired 4052 is used most of- 
ten; it’s located in the Engineering Order/ 
Processing Section. Each machine is ( v 
complemented by a 4662 Plotter. Their 
uses vary with the area that they’re in. 

The 4052 system in the Engineering Or- 
der Processing Section is used most of- 
ten; that’s where the charts are made. 
This department receives requests for 
fans to meet a certain set of specifica- 
tions — for instance, a required air flow 
volume at a given pressure and tempera- 
ture. (Other specifications, such as 
mounting and drive types, are found in 
their catalog.) A fan that meets that par- 
ticular performance intersection can be 
selected directly from the catalog tables, 
but that doesn’t tell much about the fan’s 
performance over a range of operation 
conditions. Fan performance curves are 
required to tell the whole story, curves 
that were laboriously generated by hand. 
So a program was developed to let the 
4050 System take over that task. 



Performance Curves 

The performance curves for a particular 
fan describe its operating characteristics 
over its entire operating range. While the 
full operating range is shown, emphasis 
is placed on the requested application. 
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The example in Fig. 2, for instance, 
shows the characteristics for a type 404 
fan (with a wheel diameter of 40 inches), 
moving 15,000 CFM (Cubic Feet per 
Minute). At that flow volume, empha- 
sized by the vertical line in the graph, 
the fan has a static pressure of 10 inches 
of water and requires 42.9 BHP (Brake 
Horsepower). 

Other points in the static pressure curve 
represent the total safe operating range 
for the fan. For instance, at 10,000 CFM, 
the static pressure would be about 12 
inches, while using 32 BHP. However, 
these operating areas must be balanced 
against a system requirement. Each set 
of conditions also has an operating effi- 
ciency. And mechanical efficiency trans- 
lates to energy efficiency, a topic fore- 
most in many designers minds. 

Speed-Torque Curves 



Figure 3 is a Speed-Torque Curve for a 
similar type of fan. In this graph, both 
BHP and Torque are displayed 
simultaneously, as a function of fan 
speed. This, along with wheel inertial 

(WR 2 ), is important to the designer of an Fig. 2. A Performance Curve Plot for a size 404, Series 15 GILS fan, showing the relationships of brake 
air movement system, to efficiently se- horsepower, and fan speed. 
lect the proper motor type and size for 
an installation. 

As these programs were put into use, the 
operating ease of the 4050 Series Graphic 
Computing Systems became an apparent 
asset as well. With the programs in 
place, the system can be operated by 
anyone in the group, even new staff 
members. That’s a major factor when 
you send out as many as 100 Perfor- 
mance Curve plots each day. 



Custom Shafts and Bearings 



There is other customer-directed infor- 
mation that’s generated by the 4052 in 
New York Blower’s order processing 
section. When a custom-designed fan 
system is ordered, the 4052 is put to 
work on shaft and bearing calculations. 

Using the shaft size, overhangs, and type 
and size of bearings, the program calcu- 
lates the first and second critical reso- 
nance for the fan. In addition, the sys- 
tem will calculate the required fan base. 

While these programs don’t use the 
graphic capabilities, the higher speed of 
the 4052 processor shows forth in these Fig. 3 A Speed-Torque curve for a size 491, Series 40 fan, showing the relationships of torque, brake 
calculations. horsepower, and fan speed. 
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In the Lab Too 



A separate building at New York Blower 
houses the Research Laboratory, where 
their 4051 resides, also with a 4662 Plot- 
ter. In this lab, they test all ranges of fan 
sizes and types for air flow, static pres- 
sure, mechanical efficiency and sound 
level. The system plots the results of all 
lab tests, for easy reference to test re- 
sults and quick communication of the re- 
sults to others. 



To test a fan in the lab, the fans are 
mounted in an appropriate housing to se- 
cure them for the test, and are then in- 
strumented to gather the required data. 
Small fans might be mounted in the test 
chamber wall, as shown in Fig. 4, while 
the bigger blowers might be secured in a 
free-standing fashion (Fig. 5). Of course, 
the test being performed dictates much 
of the mounting requirements. 

Plots of test data come in several forms; 
Fig. 6 is one example of a test lab plot. 
In that plot, three factors — Brake Hor- 
sepower, Static Pressure, and Mechani- 
cal Efficiency — are plotted as functions 
of air flow. Tests and test plots of this 
type help the engineering group at New 
York Blower to optimize these important 
factors in the design of new fans. 

In addition to the local processing and 
plotting with the lab’s 4051/4662 combi- 
nation, some outside timeshare process- 
ing is done as well. Finite Element Mod- 
eling is used in the design process, for 
instance, using the processing and 
software of Structural Design Research 
Corporation in Cincinnati, Ohio. To 
make use of their services, data is trans- 
ferred through the lab modem and into 
their network for analysis. The returned 
data can be plotted as well. 



Measurable Results 

Improvements are somewhat hard to 
measure in the lab, especially since they 
are using the system for things that they 
couldn’t do before. But the changes are 
easier to see in the graphs and plots the 
order processing section prepares. 

As noted earlier, graphs were formerly 
prepared through hand calculation and 
manual plotting. Using this process, each 
customer’s Performance Curve took 
about 30 minutes to produce, while the 
4052/4662 system completes them in 





/ \ 



Fig. 4. Laboratory testing of a fan used as a small 
crop dryer. 



Fig. 5. A larger, typical centrifugal type fan under- 
going analysis in the New York Blower Laboratory. 




S ^ 



Fig. 6. An example of the test plots produced in the laboratory, showing static pressure, mechanical 
efficiency, and brake horsepower, as they relate to the airflow produced. 



about 2-3 minutes. The change in pro- 
ductivity was quickly apparent, allowing 
fewer people to complete more work. 
The first system paid for itself in a year. 

In addition, there is greater acceptance 
of the validity of plots generated by the 
4052/4662 combination. No matter how 
carefully a hand-prepared plot is calcu- 
lated and drawn, it is still a hand-pre- 
pared plot, and subject to the human er- 
rors that attend such endeavors. But 
accuracy is assumed in the computer- 
drawn plots; the first hurdle is over- 
come. In the lab or in the engineering of- 



fices, the 4050 Series Graphic Systems 
continue to play an important role in in- 
creasing productivity at New York 
Blower Company. Easy use, processing 
power, and high resolution graphics are 
some of the reasons why. JSh 
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Announcing 4050D02: Statistics 
Disk-based Tests and Distributions 




Fig. 1. The 4050D02 PLOT 50 Statistics: Tests and Distributions disc-based software package. 



by Terry Davis 

TEKniques staff 

If your application requires that you 
summarize lots of numeric data, or if 
you’re a decision maker in need of statis- 
tical aid, you’ll be happy to hear about 
the newest Plot 50 software package. 
4050D02 PLOT 50 STATISTICS: TESTS 
AND DISTRIBUTIONS is a software 
tool that contains all of the functions of 
the tape-based 4050A01 PLOT 50: Statis- 
tics Vol. 1 version, and more. User in- 
terface features are enhanced. There are 
more statistical routines and more plot- 
ting options. And all of the programs are 
menu-driven and completely interactive. 

The easy-to-use nature of Plot 50 Tests 
and Distributions has been made friendly 
in several ways. The flow through the 
various levels is a natural progression 
from data entry to exploratory analysis 
to testing and analysis, or finally to plot- 



ting of results. And extensive use of the 
User-Definable Keys (UDKs) makes op- 
eration easy even for the novice. The 
UDK overlay labels the 17 keys that are 
used for data entry, editing, and storage. 

High-quality, useful plots are available 
from the analyses, including the new 
routines and options. Cumulative Histo- 
grams and Suspended Histograms are 
available, as well as Scatter Plots, Stem- 
and-Leaf Displays, and Box-and- Whisker 
plots. 

In addition, an entirely new data man- 
agement section has been added, as well 
as a new random number generation pro- 
gram. The Two-Sample F-test is also a 
new addition in this package, as well as 
the resistant line option for probability 
plots, distribution analysis, and regres- 
sion analysis. We’ll look at some of 
these in more detail, a little later in this 
article. 



Standard File Format Featured 

The PLOT 50 Standard File Format, de- 
scribed in TEKniques Vol 4. No. 1, is 
another feature; it allows data entered in- 
to, or generated by, any 4050DXX soft- 
ware package to be accessed and used 
by other “D-Series” software packages. 
Examples are Plot 50 Easy Graphing, 
and the disk version of PLOT 50 Busi- 
ness Planning and Analysis Vol. 2, both 
also described in previous TEKniques is- 
sues. The new PLOT 50 Statistics: Tests 
and Distributions can analyze data gener- 
ated by another software package, or can 
generate data to be analyzed by another 
software package. 

New Features 

The Standard File Format eases the 
interface between 4050D02 and other 
compatible software packages. But the 
most noticeable change is in the interface 
between you and the system with Statis- 
tics Test and Distributions installed. 17 
functions are defined and called through 
the User-Definable Keys, to let you se- 
lect data entry, editing, and output op- 
tions at the touch of a key, in addition to 
listing data and variables, and controlling 
program execution. The function overlay 
for the User-Definable keys is shown in 
Fig. 2, 




Fig. 2. The Statistics: Tests and Distributions 
User-Definable Key overlay. 

The entry functions allow entering data 
from the keyboard (key 1), from the disk 
(key 11), and to stop entry (key 12). The 
Editing keys allow you to Change Any 
entry, Delete Any entry, Insert Any 
entry, or Resume Entry. Output control 
consists of directing output to a plotter 
or a printer, or resetting the addresses of 
the peripheral devices. Two other keys 
allow listing variables and listing data. 
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The program execution functions on 
User-Definable Keys give you easier and 
more-direct control over program execu- 
tion. Key 6, for instance, is the Options 
function. You can press this key to list 
or perform program options. If the cur- 
rent program has no options, the follow- 
ing message will appear: 



d. 

3: 



If the current program has only one 
option, the option is executed immedi- 
ately when the key is pressed. If more 
than one option is available, the avail- 
able options will be displayed. For in- 
stance, pressing the Options key during 
the Paired t-Test will result in the 
following display, allowing you to choose 
the desired option: 



45:PL0T 58 STATISTICS: TESTS AND DISTRIBUTIONS UTILITIES 

46: 

47: 1. GET FILE 
48! 2. PUT FILE 
49! 3. DIRECTORY 
58: 4. DUPLICATE DISC 
51: 5. KILL FILE 



: 6. APPEND FILE 
: 7. LIST FILE 
: 8. CHANGE TITLE, 

: 9. COMBINE FILES 
: 10. SORT FILE 
: 1 1 . RENAME FILE 
M2. TRANSFORM DATA 



CHANGE LABELS 



63:Please choose one of the above: 



The other program control keys are the 
Repeat key, which repeats the last table 
of graph produced, and the Select Graph 
key, which is used to select the plotting 
mode. When the Select Graph key is 
pressed, one of two messages is dis- 
played: 



?: 



tt PAIRED t-TEST ** 
OPTIONS 



:l. CONFIDENCE INTERVALS 
912. PUT DIFFERENCES ON DISC 
10 : 

lllPlease choose one of the above: 1 

12 : 



The Solve key is used to rerun the pro- 
gram without recalling the menu, as 
when the variables have not changed; 
it’s faster than recalling the program 
menu. If you haven’t yet run a program 
when this key is pressed, the following 
message appears to direct you to the 
menu: 



13: 

14:S«l«ct nenu iten. Press Function Key 10 (MENU) 
15 : 



Pressing the Menu key recalls the Statis- 
tics Tests and Distributions menu, if 
you’re in the Tests and Distributions fa- 
cility: 



PLOT 58 STATISTICS: TESTS AND DISTRIBUTIONS 

** MENU »* 



22! DATA ENTRY 
23: DATA FROM KEYBOARD 

24! DATA FROM DISC 
25: 

26!DESCRIPTIUE STATISTICS 
27: 

28: EXPLORATORY PLOTTING 
29! HISTOGRAMS 
30: STEM-AND-LEAF DISPLAY 

31: BOX-AND-HHISKER PLOT 

32! SCATTER PLOT 



4i:Please choose c 



TWO SAMPLES 

T-TEST, EQUAL VARIANCE 10 

T-TEST, UNEQUAL VARIANCE 11 

T-TEST, PAIRED 12 

F-TEST 13 

SIMPLE REGRESSIONS 14 

CONTINGENCY TABLES 

2X2 15 

R X C 16 

RANDOM NUMBER GENERATION 17 



DATA MANAGEMENT UTILITIES 19 



But if you’re in the Utility facility, the 
utility menu is displayed: 



67: 

68: Sel«c t graph node is on. 



OR 



69: 

78:Select graph node is off. 
71: 

72: 



If graph mode is on, you’ll be asked to 
specify the plotting parameters; if it’s 
off, the plotting parameters will be se- 
lected by the program. 



New Utilities Too 

The statistical routines all operate under 
the Tests and Distributions facility, 
which is the default operating condition 
of the system. But a set of 13 utilities 
has been added as well, accessible by se- 
lecting program 19, Data Management 
Utilities. Some of the utility functions 
are as follows: You can GET a file from 
a PLOT 50 Standard Data Tape or 
another PLOT 50 Standard File compati- 
ble disk, or PUT a file on another Stan- 
dard Data File or another Statistics: 
Tests and Distributions disk. You can 
use DIRECTORY to obtain a listing of 
all files stored on the data disk or a 
PLOT 50 Standard Data Tape, or use 
DUPLICATE to make a copy of the 
4050D02 PLOT 50 STATISTICS: TEST 
AND DISTRIBUTIONS disk for archiv- 
ing or as another data disk. 

Other utilities let you delete data files, 
append data from one file to another, 
and list all or part of a file. You can 
change titles, change labels, combine 
variables from two files to create a new 



file, or sort an old file into ascending or 
descending order to create a new file. 
Data file names can be changed as well. 
The Transform utility lets you create a 
new file in several ways: select specific 
variables or specific observations from 
another file, or perform mathematical 
transformations on data from another 
file. 



New Programs, 

New Plots 

In addition to the easier interface to you, 
and the new utilities to ease data man- 
agement, there are several new statistical 
analysis programs in Statistics: Tests and 
Distributions. For instance, F- test has 
been added, allowing you to choose one- 
or two-sample tests. These tests display 
the sample mean, variance, degrees of 
freedom, the value of the statistic and 
the test conclusion: whether or not to ac- 
cept the null hypothesis. Let’s look at 
some of the other new programs: 

Probability Plots and Distribution Analy- 
sis is another addition to the statistics 
package. This program displays the prob- 
ability plot of any variable. It will plot 
using normal, log normal, logistic Wei- 
bull, exponential, or uniform distribu- 
tion. The fluctuations of the data points 
from the fitted line show how much the 
data deviates from the theoretical 
distribution. In addition, you have the 
option of plotting with least squares 
method or resistant line method. Figure 
3 shows a sample plot of normal proba- 
bility for water quality data. 




Fig. 3. A sample plot of normal probability for 
water quality data. 
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Regression analysis is another valuable 
program addition. Simple Regressions 
can be used to fit one of nine models to 
the data. The program will fit a staight 
line, using least squares or resistant line 
method, and plot it. Figure 4 is a sample 
of a regression analysis plot of water 
quality data, where the chosen model is 
Y = A + B * X. 




Fig. 4. A regression analysis plot of sample water 
quality data. The chosen model is Y=A+B*X. 



In addition to the new statistical analysis 
programs, the menu of plot possibilities 
is greatly expanded. This adds the flexi- 
bility of being able to explore data in a 
number of ways, graphically. Plots can 
include regular histograms, cumulative 
histrograms, and suspended histograms, 
such as the one shown in Fig. 5. In a 
suspended histogram, the distance of the 
cells from the zero line shows how much 
the sample distribution deviates from the 
normal distribution. 




Fig. 5. A sample of a suspended histogram. The 
distance of the cells from the zero line shows how 
much the sample distribution deviates from the 
normal distribution. 

Some of the other new plots include the 
Stem-and-Leaf Display, Box-and- Whisk- 
er Plots, and Scatter Plots. Stem-and- 
Leaf Displays are used to examine the 
distribution of a variable and determine 
if any values are extremely high or low. 
Similarly, Box-and- Whisker Plots give a 
quick look at the data, quickly show val- 



ues that are extreme, and therefore likely 
to be strays, out-liers, or questionable 
data points. Scatter plots are a third type 
of “check” plot, which can show the re- 
lationship between two variables, in ad- 
dition to an overall look at the data. 
Samples of these three plot types are 
shown in Figs. 6-8. 




Fig. 6. A Stem-and-Leaf display, which can be 
used to quickly examine the distribution of a 
variable. 




Fig. 7. A Box-and Whisker Plot, which can be used 
to find strays or “outliers” in the data. 
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Fig. 8. A Scatter Plot sample. Scatter Plots show 
data concentration by plotting one variable 
against another. 



More Information 

This new software package will enhance 
your system’s effectiveness, in any envi- 
ronment that requires statistical analysis 
of data. If you’d like more information 
about 4050D02 PLOT 50 Statistics: Tests 
and Distributions, contact your local 
Tektronix office. JS) 
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by D. L. Raschella 
J. R. Uchida 
J. R. Peterson 

Department of Chemistry 
University of Tennessee 
Knoxville, TN 



Background 

Our research efforts are focused on de- 
termining the basic physico-chemical 
properties of radioactive transplutonium 
elements, man-made elements produced 
in nuclear reactors. In a previous article 
we described the vital role of the 4051 
Graphic System in the acquisition and 
subsequent analysis of spectroscopic da- 
ta from these materials. This article out- 
lines the utilization of the 4051 in our so- 
lution calorimetry program. 

A calorimeter is an instrument that 
measures the quantity of heat absorbed 
or evolved in some process. Solution 
calorimetry involves determination of the 
heat absorbed or liberated when a sub- 
stance (usually a solid) dissolves in a liq- 
uid. Since transplutonium element sam- 
ples available for calorimetric study are 
typically sjubmilligram in size, our 
calorimeter was designed to detect the 
small heats involved in such dissolution 
experiments. We call our system a solu- 
tion microcalorimeter because of these 
design characteristics. 

The fundamental principles of calorime- 
try are straightforward. The experimental 
heat is measured with a carefully cali- 
brated thermometer. This calibration is 
usually done by inputting an accurately 
known quantity of electrical energy (cali- 
bration heat) into the calorimeter and 
noting the response of the thermometer. 
The thermometer’s response to an exper- 
imetal heat is then compared with its re- 
sponse to the calibration heat, thereby 
determining the experimental heat. Re- 
search calorimeters are more complex, 



depending on the type of thermometer 
employed, the precision of measurement 
desired, and the extent to which thermal 
exchange with the calorimeter’s sur- 
roundings is controlled. 

In our microcalorimeter a thermistor is 
used as the thermometer, so changes in 
temperature are detected as changes in 
resistance. A Wheatstone bridge (shown 
schematically in Fig. 1) is used to mea- 




Fig. 1 Temperature measurement system. 



sure this resistance change; the change 
of resistance is monitored by the bridge 
detector. A block diagram of the circuit 
used in the calibration of the thermistor’s 
response is given in Fig. 2. The voltages 
E h and E std are needed to determine the 
electrical energy delivered to the calo- 
rimeter during a calibration heating. 




Fig. 2 Calibration heating system. 



The microcalorimeter operates under iso- 
peribolic conditions. That is, the temper- 
ature of the microcalorimeter itself is 
allowed to change, but the temperature 
of the calorimeter’s environment is not. 
The calorimeter’s thermal excursions 
with time can be monitored on a strip- 
chart recorder or by an automatic mea- 
surement system. Our calorimetry facili- 
ty has the capability for computer-assist- 
ed data acquisition and analysis, made 
possible by the 4051 Graphic System. 

Interfacing 



Figure 3 is a block diagram of the confi- 
guration of the data acquisition/analysis 
system. The output of the bridge detec- 




Fig. 3. Block diagram of the data acquisition / 
analysis system. 



tor and E h and E std are input to a Hew- 
lett-Packard Model 3495A scanner. At 
designated times the controller (the 4051) 
commands the scanner to connect the 
appropriate channel. The 4051 then di- 
rects a Dana Model 5900 digital voltme- 
ter (DVM), through a Dana Model 55 in- 
terface, to read the voltage and send it to 
the 4051. (Dana Model 55 converts the 
voltmeter’s BCD output to the IEEE-488 
standard.) Data are collected on a real- 
time basis, with the aid of a TransEra 
Model 641 real-time clock, and are 
stored on magnetic tape using the 4924 
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Fig . 4 . Microcalorimetry facility with 405 1 system . 

Tape Unit. A 4631 Hard Copy Unit and 
a 4662 Interactive Digital Plotter provide 
permanent records of the output. All in- 
terfacing is achieved through the GPIB. 
Dr. Raschella is shown with the 4051 and 
microcalorimeter systems in Fig. 4. 

Data Acquisition and Analysis Soft- 
ware 

The program is divided into two parts: 
acquisition, and analysis/evaluation. The 
input section is stored in one program 
file; the analysis/evaluation section is lo- 
cated in several files, with overlaying 
achieved by depressing appropriate 
User-Definable Keys (UDKs). 

Data Acquisition. The overlay card for 
the data acquisition routine is shown in 
Fig. 5. The program tape is inserted into 
the 4924 tape drive, while the data tape 
is placed into the 4051 tape drive. New 
data tapes are automatically MARKed 
by pressing UDK 15. After entering the 
sample identification number and date, 
the data acquisition process is begun by 
depressing UDK 1. The scanner then 
samples the bridge detector output at a 
regular rate. Reading rates may be 
altered from about three readings per 



e digital calorimetry™ 6 ' input fil6 ' 




Fig. 5. User- definable keys’ overlay^card for data 

acquisition. 



second to one reading every three sec- 
onds by pressing UDKs 6 through 10. 

The collected experimental data consist 
of bridge calibrations (commonly re- 
ferred to as chart sensitivities), calibra- 
tion heatings, and the dissolution 
reaction heat of the sample. Typically, 
several chart sensitivities are deter- 
mined, interspersed with calibration 
heatings, before the dissolution reaction 
heat of the sample is measured. Then 
several more chart sensitivity and cali- 
bration heating determinations are made 
following dissolution of the sample. (For 
further information regarding calorimetry 
and the details of the experimental pro- 
cedures, refer to reference 2.) Chart sen- 
sitivities provide the relationship 
between changes in voltage (AV) and 
changes in resistance (AR) and are a 
measure of the combined sensitivites of 
the thermistor and bridge detector. They 
are performed by changing the decade 
resistance setting (see Fig. 1) and noting 
the resultant change in bridge detector 
output. Decade resistance values are in- 
put from the 4051 keyboard after de- 
pressing UDK 3. Calibration heatings ef- 
fectively calibrate the thermistor in units 
of heat (joules) per decade resistance 
change caused by the heat input. 

The 4051 is prepared for a calibration 
heating by depressing UDK 5, at which 
point E h and E st( i, as well as the bridge 
detector channel, are scanned sequential- 
ly at a rapid rate. The actual input of 
heat as electrical energy occurs in a 
fixed time period and is initiated by a 
switch on the calorimeter. After the cali- 
bration heating the values of E h 
[EMF(l)] and E sW [EMF(2)] are dis- 
played, and the new bridge resistance 



setting must be entered. All events must 
be identified: I for initial entry, C for 
chart sensitivity, H for calibration heat- 
ing, and R for the dissolution reaction. 
Examples of the 4051 ’s CRT displays 
during an experiment are given in Fig- 
ures 6a, 6b, and 6c. User-definable key 
1 1 temporarily stops data collection but 
not the clock, while UDK 12 stops both 
the data collection and clock and closes 
the data file. User-definable key 12 is 
useful if an interruption in data collection 
is prolonged, such as for the adjustment 
of some equipment. The end of the ex- 
periment is designated by UDK 13, 
which closes the current data file and 
stores the sample identification number 
and a summary of information about the 
chart sensitivites and heat inputs in file 1 
of the data tape. Analysis of the data is 
begun by pressing UDK 16. 



Enter Sanple 


ID: Cn<248> 


207-6-E2 




Enter Pate: ? 


-25-78 






Ready. Press 


User-Definable 41. 




TIME-103.2 


R-426419 






Chart, Heat 


, Reaction, 


or Initialize: I 




TIME-521.5 


R-426355 






Chart, Heat 


Reaction, 


or Initialize: C 




£HF<I>: 








1.9373 


0 


0 


0 


0 


0 


0 


0 


EMF<2>: 








1 . 945598 


e 


0 


0 


0 


8 


0 


0 


TIME -978.6 R-426218 






Chart. Heat, 


Reaction, 


or Initialize: H 




TIRE-1331. 6 


R-426158 






Chart, Heat, 


Reaction, 


or Initialize: C 




EHF<1>: 








1 . 9373 


1.9373875 0 


0 


0 


0 


0 


0 


EHF<2>: 








1.943598 


1.9456 


0 


0 


0 


0 


0 


0 


TIHE- 1736. 9 R-42601? 






Chart, Heat. 


Reaction, or Initialize: H 




TIME-2419. 3 


1-425957 






Chart, Heat, 


Reaction, or Initialize: C 




£HF<1>: 








1 . .9373 


1.9373875 1.9373 


0 


0 


0 


0 


0 


EMF C2>: 








1 . 945398 


1.9456 


1 . 945596 


> 0 


0 


0 


0 


0 


TIME-321 1.0 


R-425816 






Chart, Heat 


, Reaction, 


or Initialize: H 




TIME-4689.3 


R-425756 






Chart, Heat 


, Reaction, 


or Initialize: C 




TIME-6258.6 


R-425606 






Chart, Heat 


Reaction, 


or Initialize: R 




Fig. 6., a, b, c 


Sample CRT displays during data 



acquisition. 
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Data Analysis/Evaluation The overlay 
card for the data analysis/evaluation rou- 
tine is shown in Fig. 7. User-definable 
key 5 is used to plot the voltage- time da- 
ta. An example of such a plot is given in 
Fig. 8, where chart sensitivities occur 
around 450 and 1180 seconds and cali- 
bration heatings at about 750 and 1600 
seconds. User-definable key 1 initiates 
the chart sensitivity analysis routine. 
This extrapolates the drift line segments 
prior to and following the resistance 
change towards an “imaginary” line that 
intercepts the voltage-time curve halfway 
through the event which caused the 
change in resistance. Then the distance 
(AV) between the points at which the ex- 
trapolated drift line segments intercept 
the “imaginary” line is determined. The 
drift line segments to be extrapolated are 
selected by the operator by using UDKs 
6 through 10. This procedure is illustrat- 
ed in Figs. 9a and 9b. 




Fig. 7. User-definable keys’ overlay card for data 
analysis /evaluation . 




Fig. 8. Voltage-time plot of some typical calori- 
metric data. 




1751 "tT 3 e"~iM 8 iTse'iiSe 'lire ifiea tBa "1200 '“1210 '1220 1230 

r CHART CAl.IBRATIGM 

Identify Line Segments to Extrapolate R-Change: 




Oil Terence « -0.8933391396879 



Fig. 9., a, b Analysis of chart sensitivities . 



Calibration heatings are evaluated in a 
similar manner. Drift line segments be- 
fore and after the heat input are extrapo- 
lated towards a point midway into the 
heating period. The “ imaginary” line 
for calibration heatings and the dissolu- 
tion reaction heat is known as “T M .” 
(An example of its determination is given 
later.) The distance between the extrapo- 
lated drift line segments at T M is then 
measured. User-definable keys 3 and 6 
through 10 are used in this case. The ex- 
trapolations are illustrated in Figs. 10a 
and 10b. The beginning (B), end (E), and 
T m are automatically marked by the pro- 
gram. 




Fig. 10., a, b Analysis of calibration heating. 



A stepwise illustration of a determination 
of T m is given in Figs. 11a through lid. 
Using UDKs 2 and 6 through 10, line 
segments are identified (Fig. 11a). These 
line segments are extrapolated inward, 
tangentially from the inner point chosen, 
in discrete steps, and the areas between 
the extrapolated lines and the experimen- 
tal curve are calculated for each step, as 
shown in Figs, lib, 11c, and lid. These 
areas are determined by a modified New- 
ton’s integration technique. (Alternative- 
ly, a nonlinear extrapolation of the 
chosen line segment might be used.) T M 
is the one point at which the areas be- 
tween each extrapolated line, the experi- 
mental curve, and the T M line are equal. 
Once T M is found, the drift line segments 
are extrapolated to T M as before, and the 



difference between the segments at T M is 
determined (see Figs. 12a and 12b). 
Upon concluding the data analysis, de- 
pressing UDK 19 yields a summary of all 
pertinent information. This information is 
used in subsequent calculations to deter- 
mine the heat of solution of the sample. 



3i'<S 364ir~ 3636" 3663 5670 

* Th DETERMINATION 

Identify Lin* Segxents to Extrepol 

Enter StS traight LIn«> or CtCurwed L 


• ■3677 'S6B4 ~S691 "5698 STOT 

ate Heat Effect: 

Inc) Extrapolation: S 


■ 


3677 3684 5691 3698 5703 

ate Heat Effect: 

ine) Extrapolation: S 


/ 

3642 3649 5636 3663 3670 5677 5684 5691 3693 3, ■’03 

X Tn DETERMINATION 

Identify Line Segments to Extrapolate Heat Effect! 

Enter s<strat»ht Line* or CCCurved Line) Extrapolation". S 


u 

5643 5649 5636 5663 567 

* Tn DETERMINATION 

Identify Line Segnent* to Ex traps 

areal * 0.285737733281 
area2 - 8. 283383923906 
Tn » 3672.3898623 

Enter SiStraight Line) or CtCurved 


r 

0 3677 5634 5691 5698 5703 

late Heat Effect: 

Line) Extrapolation: 3 



Fig. 11., a, b, c, d Determination off M. 
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* DR !P1 LIHE EXTRAPOLATION 

Idantify Lin« 3«»rt«rit* ?o Ext'-apelat* Heating 
Correction Factor - 8. eel 18383SM97? 

Reaction or Heat? R 



Fig. 12., a, b Analysis of dissolution reaction heat. 



Conclusions 



An appreciable number of extrapolations 
and integrations are required in evaluat- 
ing calorimetric data; the 4051 is 



particularly helpful in these aspects of 
data analysis. Not only can the 4051 re- 
duce the time required for evaluation of 
the experimental data, but it can also 
eliminate “eye-balling” errors inherent 
in “manual” extrapolations and integra- 
tions. The 4051 has also increased our 
flexibility, by providing the ability to 
perform the necessary extrapolations on 
a nonlinear basis, as well as on a linear 
basis. The 4051 offers unsurpassed accu- 
racy and ease in this regard; its capabili- 
ty for regression analysis and integra- 
tions is invaluable. Jgp 
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INPUT / 
OUTPUT 








Jerry W. Anderson, Phillips Chemical 
Company of Bartlesville, Okla, has several 
questions. Howard Sanders, Technical 
Support Specialist at Tektronix replies. 

Plotter/Printer 
On or Off 

Is it possible to determine, under 4052 
program control, whether or not the plot- 
ter and/or printer is turned on? (The mo- 
tive here is to output a plot to the plotter 
if it is turned on; otherwise output it to 
the screen.) 

The 4052 by itself cannot determine if a 
plotter or printer is turned on before it 
attempts to output (or input). If the out- 
put is to the plotter and it isn’t on, you’ll 
receive an error message and your program 
will halt. If the output is to the printer 
which isn’t on, the 4050 system 
will “go busy” waiting to output. How- 
ever, a new 4050 ROM Pack, the 4052R09 
Real Time Clock ROM Pack* provides the 
facility to interrupt or continue a program 
after a prescribed amount of time has elaps- 
ed. Taking advantage of this ability, the 
following routine recovers your program if 
the plotter isn’t on. 




The “ONTIME” command tells the ROM 
Pack to branch to line 84 after a specified 
amount of time. If the program is running 
when the time is up, it branches and when 
it encounters a RETURN, returns to the 
statement following the one it executed prior 
to branching. If a program isn’t running, 
the ROM Pack starts the program execut- 
ing at statement 84 (a RETURN in this case 
would return it to idle mode). 

Applying this, the above routine branches 
to line 84 after .21 seconds. The MOVE 
in statement 130 will generate an error if the 
plotter isn’t on and the program will halt; 
the variable, F, remains 0. If the plotter is 
on, F is set to 1 and the program continues. 
Now, when the ROM Pack instructs the 4050 
system to branch, the value in F at line 84 
tells whether the program was executing. 
If it is 1, it simply returns to the program. 
If it’s 0, the device address is set to 32 and 
the program directed to the plotting routine. 



However, there is no way with or without 
the ROM Pack to recover from the busy 
state, and determine if the printer, whether 
interfaced through Option 10 or Option 1, 
is turned off. 

Enabling BREAK on Option 1 

Using the Data Communications Inter- 
face, is it possible to enable the “BREAK” 
key without also enabling the “PAGE- 
FULL” break? 

This is accomplished by disabling the 
PAGEFULL break. To do this set the first 
parameter in CALL “MARGIN” to 0. 

From Bill Fowee at Proctor & Gamble in 
Cincinnati comes the next question; again 
Howard Sanders responds. 

DIR @1:0, “@” 



niques. I found last month’s “TLI @1:” 
useful. How do you DIR @1:0, “@”? 

The 4050 Series can’t directly output the 
directory of the 4907 to the GPIB, however, 
it can output to the tape. Thus, using the 
tape as an intermediate step and a small 
program to input and output the directory 
tape, the directory can easily be sent to 
the plotter. 



100 REM Routine to put di*k directory on 
110 REM tape and then to plotter 
120 INIT 
130 FIND 1 

148 DIRECTORY @33:0,"*" 

150 CLOSE 
160 FIND 1 

179 ON EOF <0> THEN 210 

180 INPUT §33:8$ 

190 PRINT §l:A$ 

208 GO TO 180 
210 END 



USMAP/SEG67 

USMAP/SEG76 

USMAP/SEG85 

USMAP/SEG94 

USMAP/SEG5 

USMAP/SEG 1 4 

USMAP/SEG23 

USMAP/SEG32 

USMAP/SEG39 

USMAP/SEG4 1 

USMAP/SEG48 



The 4907 directory can be output directly 
to the Option 10 Interface (RS-232-C) and 
thus to the printer by simply typing: 




*For use only with the 4052 and 4054 systems. Thanks for your new column in TEK- 
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The CALL 6 ‘FILE” command may also be 
used. A full file status will be generated 
for each file in the library designated and 
then stored in a string variable. The 
variable may then be output over the GPIB 
to the plotter. 




The limitation, however, is plotter paper 
area. If you ask for a total file directory 
dump, e.g., CALL “FILE”, 0, “<§>”, A$, 
one sheet of plotter paper may not be 
enough to hold it all and you can’t change 
paper without losing some of the string. 



Fitting Data 

Dan Taylor, Design Engineer at Tektro- 
nix and TEKniques’ Technical Editor 
writes: Twice recently a question has 
come to me about data fitting. The cus- 
tomers had a set of (X^YO observations 
and a model Y = F(X,C) and wanted to 
find the “best” C to fit the data (where 
C is defined as equal to one or more 
“constants” which are chosen to fit the 
data). 

There are several different programs in 
PLOT 50 Statistics Vol. 4 that will do this 
(find C). Marquardt’s method was used 
by both of these customers. However, what 
they really wanted was Y = F(X,C) with 
the additional constraint that C is greater 
than or equal to 0. 

This is quite possible and easily done; just 
code the model in BASIC as Y=F(X,ABS 
(C)). For example, the model Y=aE bx would 
be coded as: 

1180 Y*ABS<C< 1 ) )*EXP<ABS<C<2) )*X< I > > 



the Technical Support Specialists. Their 
answers may save you a telephone call 
to your Systems Analyst. 

ON SRQ and POLLing 

When using a 4050 with an LSI (Digital 
Equipment), the second SRQ doesn’t get 
trapped coming from the LSI although a 
POLL has been executed. 

The ON SRQ THEN . . . statement traps 
the first SRQ, but to arm the 4050 System 
to respond to another SRQ, a RETURN 
must be incorporated somewhere after your 
POLL routine to complete the SRQ loop. 
When the RETURN is executed, the ON 
SRQ THEN . . . will respond to the second 
interrupt, and so on. For a thorough dis- 
cussion on the handling of interrupts from 
peripheral devices on the GPIB, see “The 
POLL Statement” Section of your 4050 
Series Graphic Systems Reference Manual. 



110 ON SRQ THEN 1970 
120 REN Begin Main program 



1970 POLL M,WJ5|8S3 

1980 GO TO M OF 2000,3080,4000 

I 

2000 REN Service Routine for Device Number 5 



2990 RETURN 

I 

I 

3000 REM Service Routine for Device Number 8 

1 

I 

3990 RETURN 

I 

I 

4088 REN Service Routine for Device Number 3 



4990 RETURN 



Tape Error When OLDing 

I saved a program on tape in ASCII for- 
mat. When I attempted to OLD the pro- 
gram back, a tape error occurred in the 
middle of the file. How can the rest of 
the program be retrived? 



Error Message MAG TAPE READ ERROR - MESSAGE NUMBER 53 
indicating bad spot 



Key in INPUT @33: A# 

See what's in A$ 

A* Data Bases, enter! FILEINIT 

Do it again to INPUT <?33! A# 
get line number A$ 

320 PRINT " tor Mapping Demonstration, 
enter: SEGUSMAP 



Key in statement 320 PRINT " for Mapping Demonstration, 
enter: SEGUSMAP" 

Key in next line # 330 REM 
Append rest of APPEND 330 
program 



| Data Bases, enter : FILEINIT" 320 : : : 

t 

File Pointer 

(File pointer will be placed after bad spot 



4907 AUTOLOAD 



3 INIT 
3 PAGE 

5 C*«“4-JUL-76" 

5 D*«"tl:li:ii" 

3 PRINT " 

3 PRINT ’ 

3 PRINT -ENTER TODAY'S DATA (format: DD-Mon-YY, «.g.: 4-JUL-76): “s 
3 INPUT A* 

5 IF A»<>“" THEN 160 

I A*»C* 

3 PRINT "ENTER CURRENT TIME (format! HH:MM:SS, «.g.: 18:58:46): 

3 INPUT B* 

5 IF B*<>"» THEN 180 
S Bt»D* 

3 ■ 

3 A»-A*tB* 

3 CALL "SETT IN", A* 

3 CALL "MOUNT ",0,"" 

3 PAGE 

3 PRINT ’ for Napping Demonstration, enter: SEGUSMAP" 

3 PRINT "I for Data Sorting, enter: SORTING.!" 

3 INPUT A* 

3 CALL "USERLIB", "DEMOS" 

3 OLD A* 



Salvaging Accidentally Re-MARked Files 

How do I recover files which have become 
inaccessible because the tape was in- 
advertently re-MARked? 

Two programming tips for recovering such 
files on a 4051 have been published in TEK- 
niques Vol. 2, No. 3 and Vol. 3, No. 3 
(pages 39 and 75 in the new Programming 
Tip Handbook) and one tip for recovering 
these files on a 4052/4054 in TEKniques 
Vol. 3, No. 8 (page 102 in the Programming 
Tip Handbook). Jg ft 



instead of 

1100 Y«C<1>*EXP<C<2>*X<I>> 

This will work for any sort of program 
which fits (X,Y) data through a user cod- 
ed Y = F(X,C) function. 

Howard Sanders has compiled three of 
the most frequently asked questions of 



OLD in your file. When the error message 
is received, key in INPUT A$. You will 
probably get the last part of a statement. 
So do another INPUT A$ to get the next 
line number that’s retrievable. Now key in 
that statement. Then type in the next line 
number followed by REM; APPend the rest 
of the program to that line. This procedure 
can be repeated if more errors occur. 
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Editor’s 

Note 




New Catalog Arrives 

With your last TEKNIQUES issue, you 
should have found a copy of the new 
4050 Series Applications Library catalog. 
The catalog has added 67 new programs 
to the library since last year’s catalog. 
Look through it; you’ll probably find 
some you can put to work in your appli- 
cation. 

Additional copies are available for the 
asking, from the Applications Library of- 
fice serving you. Their addresses are lo- 
cated at the back of each issue. The Ap- 
plications Library staff thanks all of you 
who contributed programs to the library, 
for making these programs available to 
others. 

The Programming Tip 
Handbook is Here Too! 

Programming Tips have proven to be one 
of our most-read features in TEKniques 
over the past few years. They’ve proven 
to be such a valuable reference for most 
of our readers, in fact, that we’ve com- 
piled all of the Programming Tips and 
BASIC Bits from Volumes 1, 2, and 3, 
and put them together in a Programming 
Tip Handbook. 

The Programming Tip Handbook is avail- 
able through the Applications Library of- 
fice serving you. It’s listed in the Re- 
source Materials section of the new 
catalog; its Abstract Number is 51/00- 
7004/0. And again, the Applications Li- 
brary staff would like to thank those of 
you who contributed your valuable Pro- 
gramming Tips. They’ve proven very 
helpful to many users, and made our 
handbook possible. 



Back Issues are Gone 

TEKniques has made a point of keeping 
some back issues available for those of 
you who may have lost one, or want to 
pass a copy on to someone else. How- 
ever, we will no long provide back issues 
of copies from Volumes 1 through 3 
(1977-1979). But the information is still 
available. Many back issue requests 
were from those who wanted all of the 
Programming Tips and Basic Bits; these 
are all contained in the Programming Tip 
Handbook. Others want a copy of an ar- 
ticle relating to a particular application; 
these are now available in article reprint 
collections from the first three* volumes 
of TEKniques. (See the following note.) 
And, of course, all of the Applications 
Abstracts are found in the new Applica- 
tions Library Catalog. 



Application Article Reprints 

To continue to meet the needs of users 
who may need to review an application 
article from one of the previous issues of 
TEKniques, we’ve made some reprints 
available. We collected the application 
articles from our first three volumes, and 
put them together by application area. 
We currently have five sets of article re- 
prints available, in the following areas: 



Engineering and Design AX-4449 

Mapping AX-4460 

Data Acquisition and Analysis . .AX-4450 
Business Graphing 

and Reporting AX-4451 

Peripherals and ROM Packs . . . .AX-4452 



If you have need of an article from an 
issue in Vol. 1, 2, or 3, contact your lo- 
cal Tektronix office, or the Applications 
Library office serving you. They can 
provide the reprint volume that contains 
the article you want. 
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New Contest Coming 

Response to the last Applications Li- 
brary contest was tremendous. The en- 
tries provided several new and useful 
programs for the library. So we’re going 
to have another contest very soon. 
Watch the next issues of TEKniques for 
details. 

Wanted: 

Your Programs 

We’re always looking for more programs 
to add to the 4050 Series Applications 
Library. This helps us keep the library 
growing to better serve more of you, and 
provides a medium for you to share your 
programs with others who may be able 
to use them too. And you get three free 
programs for each program you send in 
to the library, so everybody wins. 

Programming Tips Too 

The Applications Library staff would al- 
so like to remind you that we’re looking 
for those valuable Programming Tips and 
BASIC Bits. You’ll receive any one of 
12 programs from the library in exchange 
for your tip. Details are listed at the end 
of the BASIC Bits column in each issue. 

JED 
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PROGRAMMING 

TIP8 



Checking Numeric Strings 
for Valid Numbers 



by Jim Gish 

Tektronix, Inc. 

Irvine, CA 

There are times when you may want to 
input numbers from the keyboard into a 
character string. One instance might be 
an indefinite loop to input numbers from 
the keyboard until the operator responds 
with “END” instead of number. How- 
ever, if the program attempts to convert 
“END” into a number, it will abort the 
program. To avoid this, use a numeric 
parser. A parser, in this sense, deter- 
mines if an input string contains a valid 
number without actually converting it to 
a numeric quantity . 

If the parser determines the string has a 
valid number in it, then it’s safe for the 
program to convert the string into a 
number. If not, the program can ask the 
operator to try again or check to see if 
the operator wants to exit. 

This particular parser also edits the char- 
acter string, removing any imbedded 
blanks or commas. It will accept integer, 
real or scientific notation. 



The string to be parsed is 0$. 0$ will 
have a length equal to the value stored in 
variable 13. The length of 0$ must never be 
greater than 13. When the subroutine is 
finished a numeric constant will be right 
justified in 0$ with blank fill from X$. If 
the parser doesn’t find any numeric char- 
acters KO will be set to one as a warning 
to the calling program. 



1000 REM SUBROUTINE TO UERIFY NUMERIC STRING 
1010 REM STRING PASSED IN 0* 

1020 RESTORE 1030 
1030 DATA 010,0,0.0,0 
1040 READ K0,H4,K1,K2,K3,K4 
1050 DELETE J*,U* 

1060 H1>LEN(0$> 

1070 IF H1»0 THEN 1490 
1080 DIM J$(H1),U$<1) 

1090 J*«*" 

1100 U*=“" 

1118 FOR H2»l TO HI 
1120 Y*«SEG<0*,H2,1> 

1130 K4«P0S< "0123456789. +-E, ",Y*,1> 

1140 IF K4<-0 THEN 1498 

1150 REM CASE 1 - Y* ■ A COMMA OR BLANK "IGNORE IT- 
1160 GO TO K4-10 OF 1440,1400,1480,1330,1230,1230 
1170 REM Y* IS A NUMERIC CHARACTER 
1180 K3-1 

1190 REM ADD CHARACTER TO OUTPUT STRING <J»> AND INC COUNTER <H4> 
1200 H4-H4+1 

1210 J*»REPm,LEN<J*> + l,0> 

1220 U*»Y* 

1230 NEXT H2 

1240 REM CHECK J* FOR ZERO LENGTH OR NO DIGITS 

1250 IF K3-0 OR LEN<J*>O0 THEN 1490 

1260 REM PACK NUMERIC STRING RIGHT JUSTIFIED 

1270 REM X* CONTAINS A STRING OF BLANKS AT LEAST 13 LONG 

1280 0*«J* 

1290 H*»SEG<X*,1,I3-LEN<0*>> 

1300 0*-REP<W$, 1,0) 

1310 RETURN 

1320 REM CASE 2 - Y*« "E" 

1330 IF LEN<U*>»0 THEN 1490 

1340 IF ASCCU»)-43 OR ASC<U*>-45 THEN 1490 

1350 Kl-l+Kl 

1360 IF K 1 > 1 THEN 1490 

1370 K3=0 

1380 GO TO 1190 

1390 REM CASE 3 - Y# = A PLUS OR MINUS SIGN 

1400 IF U*- ,,B THEN 1190 

1410 IF ASC<U*><>69 THEN 1490 

1420 GO TO 1190 

1430 REM CASE 4 - Y* » A DECIMAL POINT 

1440 K2=K2+1 

1450 IF K2>1 THEN 1490 

1460 IF Kl->1 THEN 1490 

1470 GO TO 1190 

1480 REM *** ILLEGAL ENTRY OR END OF INPUT 
1490 K0“1 
1500 RETURN 



Ed. Note: Thisroutineusesmany 
scratch variables. An alternative method 
is to concatenate the input string with an 
absurd number and then test the value of 
0$, e.g., 



100 REM 0**,", 3. 6492871 41697E304" 

lie REM K0«<UAL<O$)< >3. 64928714 1697E+304) 



Where K0= 1 means okay and K0=0 
means error. This is similar to the pro- 
gramming tip in Vol 4. No. 1 ‘Use VAL 
Function to Avoid Errors. 4 



APPEND to Recover 
from Alternate Delimiters 



by Dan Taylor 
Tektronix, Inc. 

Wilsonville, OR 

TEKniques Vol. 2 No. 7 and Vol. 3 No. 
1 (pages 56 and 63 in the Programming 
~~ Tip Handbook) contain programming tips 
to recover ASCII programs SAVED with 
alternate delimiters, namely CR/LF. 



Although you can’t OLD in a program 
containing CR/LF (PRI @37,26:1) as its 
delimiters, you can APPEND it if your 
environmental parameters are set to 
CR/LF. Input the first line to get the line 
number, type a REM statement for that 
line number, then APPEND in the pro- 
gram. 



PRI937, 26: 1 
FIND 1 
IMP 033: A* 

A* 

100 INIT 

180 REM 
FIND 1 
APPEND 100 

PRI 037,26:0 
LIS 

100 INIT 

110 PRINT -TEST PROGRAM* 
120 END 
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Input Without Echo The secon d routine is used for single key 

response. The examples which call it ex- 
by S. Schicktanz ecute a GOTO or GOSUB branch on the 

Physics Department response. 

University of Munich 
Germany 

The following subroutines allow a user to 
key in a response without echoing the re- 
sponse to the screen. The first routine is 
used for multiple key entries for numbers 
or strings. Two examples for calling the 
subroutine are shown. 



Note: Pressing the RETURN key 
initiates the branch to the first line num- 
ber following OF-in this case, line 400. 

Editor’s Note: Mr. Schicktanz has touch- 8, 10; Vol.2, Nos. 3, 6; (pages 16, 18, 24, 
ed on areas of several previous Program- 43, 92 in the Programming Tip Handbook); 
* T$ cannot hold more than 72 characters ming Tips— no echo, branching, checking and Vol. 4, Nos. 1,3; this issue also con- 
unless redimensioned, and it must con- for valid input and GOSUB/GOTO. For tains a tip which checks for valid numeric 
tain at least one character upon return. a review, see TEKniques Vol. 1, Nos. 7, input. 

Determination of a 
Graph’s Upper Bound 

by Leonard Weitman 
Tektronix, Inc. 

Beaverton, Oregon 

While programming a 4050 series Graph- 
ics System to display data via line and 
bar graphs, a consistent problem arises: 

What will be the largest label on the ver- 
tical axis (the top of the graph)? Often, 
data from one graph to the next will vary 
by several orders of magnitude. For this 
reason, a general utility routine for deter- 
mining the maximum graph value must 
be able to handle a large variety of mag- 
nitudes of data. Not only should the 
maximum value on the vertical axis suit 
the data’s order of magnitude, but it 
should also be divisible by some “nice” 
number so that the graph will have a 
pleasing appearance. The following sub- 
routine calculates the maximum graph 
value and the number of major horizon- 
tal tic marks that will appear on the 
vertical axis. The routine requires that: 

1) The largest data value has been de- 
termined (“L”). 

2) The lower bound of the graph is ze- 
ro. 




The number of major horizontal tic 

marks on the vertical axis will always be i 

3, 4 or 5. The increment of each of these 
tic marks will be evenly divisible by 

some pleasing quantity, such as multiples , 

of 5 or 10. Also, the maximum data point 
will appear in at least the upper third of 
the Graphics System display window. 

Although this routine assumes a lower 
bound of zero, similar logic can be used 
to determine both a variable upper AND 
lower bound. 
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Multiple Dimensioned 
Arrays 

by Nathan Oxhandler 
Tektronix, Ine. 

Santa Clara, C A 

and Dan Taylor 
Tektronix, Inc. 

Wilson ville, OR 

There are times you need to work with 
data arrays of three or more dimensions. 
While the 4050 Desktop Computers pro- 
vide for two-dimensional arrays, it’s very 
easy to develop routines that allow you 
to simulate arrays of almost unlimited di- 
mensions. 



Notice that to simulate DIM A 
(N1,N2,N3) and access it in the manner 
A(I,J,K) . . ., you need merely 

DIM A<N1,N2*N3> 

(statement 150) 

and access it as 

A<I, <J-i>*N3+K> = 

(statements 180-240) 

Multi-Dimensioned Array in Two Di- 
mensions 

Based on the foregoing, let’s extend the 
number of dimensions. 




Notice the subscript in statements 210 
and 300 which allows you to access your 
arrays in the normal I, J, K, loop index- 
ing. 



Memory 

One of the first things to look at is the 
amount of memory numeric data takes. 
Each single data item takes up eight 
bytes of memory. A single-dimensioned 
array takes eight times the number of 
elements in the array plus 18. Example: 
DIM A(20) requires 8 X 20 + 18 bytes. 
A two-dimensional array takes eight 
times the number of elements in the first 
dimension, times the number of elements 
in the second dimension, plus 18 bytes. 
Example: DIM A(20,30) requires 8 X 20 
X 30 + 18 = 4,818 bytes. As you can 
see, no matter what type of array you’re 
working with, you’ll need to be con- 
scious of the amount of memory it re- 
quires. 



Three-Dimensional Data Array in 
Two Dimensions 

The following routine simulates a three- 
dimensional data array using a two-di- 
mensional array. The application is sim- 
ple, but demonstrates the subscript han- 
dling to achieve the simulation. 




100 INIT 

110 REM MULTI-DIM ARRAY IN 2D 
120 REM X-NUMBER OF DIMENSIONS 
130 X=4 

140 REM H CONTAINS THE LENGTH OF THE DIMENSIONS 
150 DIM NCX> 

160 DATA 4, 3, 2, 5 
170 READ N 

180 REM D IS LENGTH OF 2ND DIM 
130 D=1 

200 FOR 1=2 TO X 
210 D=N(I>*D ‘ 

220 NEXT I 

230 DIM ACNC1>,D> 

240 A=0 

250 R=RND<-0. 5) 

260 FOR 1=1 TO NCI) 

270 FOR J=1 TO N<2> 

280 FOR K=1 TO NC3) 

290 FOR L=1 TO N<4> 

300 ACI, CCJ-1)*NC3)+CIC-1>>*NC4)+L>»INTCCRNDC1)+0.005>*100> 

310 NEXT L 
320 NEXT K 
330 NEXT J 
340 NEXT I 
350 PRINT A 

360 REM TO ACCESS ALL N3 ELEMENTS IN THE N2 DIMENSION IN N1C3) 
370 FOR 1=3 TO 3 
380 FOR J=1 TO NC2) 

390 FOR K=1 TO NC3) 

400 FOR L=1 TO N<4> 

410 PRINT ACI, CCJ-1)*NC3) + CIC-1>>*NC4>+L> 

420 NEXT L 
430 NEXT K 
440 NEXT J 
450 NEXT I 
460 END 



To DIM A(N1,N2,N3,N4) and access it 
A(I,J,K,L) = • • you simply 

DIM A<N1,N2*N3*N4> 

(statements 180 — 230) 



and access it by extending the subscript 



A<I,<<J-l)*N3+<K-imN4+L> = . . 

(statements 260 — 340) 



For five dimensions: 



DIM ACN1,N2*N3*N4*N5> 



and access it: 



A(I, <<<J-l>*N3+<K-imN4+<L-l>>*N5+M> 



and so on. 



Multi-Dimensional Array in One Di- 
mension 



A three-dimensional array simulated in 
one dimension is illustrated in the follow- 
ing routine: 



The subscript could easily be incorporat- 
ed as a record number in a disk file. 



As the dimensions increase, it’s a simple 
matter to expand the array subscript to 
reflect them: 




Of course, the amount of data you may 
have in your 4050 computer is limited by 
the size of your memory. However, it is 
possible to simulate very large arrays of 
many dimensions using the 4050 Series 
and the 4907 File Manager. This will be 
covered in a subsequent issue of TEK- 
niques. JSjh 
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String Concatenation 

by Stan Jensen 
State of California 
Air Resources Board 
Sacramento, CA 

The following program and sample run 

demonstrate an alogrithm for string con- “ 

catenation which is a major improvement Editor’s Note: This algorithm is intended 
over the Programming Tip on the same for building a string of the same charac- 
subject in TEKniques, Vol. 1 No. 1. ter. The length of C$ must be equal to 1. 



3 REM 
3 INIT 

3 DIM B4<18)»C$(l)fE$(i8) 

3 C*="*" 

3 L=2tl4+1 

3 PRINT "iJJJJJIILength - "JL 
3 PRINT "JJJJJIIFAST METHOD 
3 CALL "TIME'SB* 

3 DIM S*<L> 

3 S$="" 

N=INT(LOGa)/LOG<23> 

S$=C$ 

IF L=1 THEN 310 
FOR 1=1 TO N 
S$=REP(S$j LEN(S$>+1 )0> 

NEXT I 

IF LEH<S4>=L THEN 310 
DIM T4(L-LEN<S$:0 
T*=SEGCS$,1,L-LEN<S4)> 
S$=REP(T4,LEN<S4)+1,0> 
DELETE T4 
CALL "TIME" j E4 
PRINT "USTART G ";B4 
PRINT "USTOP G 11 i E$ 

PRINT " JJJJJIISLOW METHOD 
CALL " TIME", B4 
DIM S4<L> 

84="“ 

FOR 1=1 TO L 
84=REP(C$i LENCS4 > + l > 03 
NEXT I 



Length - 16365 



FAST METHOD i 

START G 20-HOU-79 11:14:23 
STOP G 20-NOU-79 11:14:25 



SLOW METHOD : 

START G 20-NOU-79 11:14:25 
STOP G 20-NOU-79 11:18:34 



RUN 



Bar Shading 

by John J. Brophy 
Pentagon 
Washington, D.C. 

The following subroutine quickly draws 
and shades bars. The driver program 
(statements 100 — 280) would have to be 
optimized for a user’s requirements. In 
this case up to four bars are provided for 
varying in height. L, R, B, T are coordi- 
nates of the LEFT side, RIGHT side, 
BOTTOM and TOP, respectively. The 
stripe code, P, is defined as follows: 

P > 0 P vertical lines/bar (not counting 
sides) 

P = 0 No shading 
P = — 1 Lines slanted to upper left 
P = —2 Lines slanted to upper right 
P > 0 P horizontal lines/inch 




Programming 
Tip Exchange 

Send in your programming tip. Any one of 
the following 4050 Series Applications Li- 
brary programs* will be yours when if s 
published. Simply jot down a brief descrip- 



tion of the function , the code, and your 
choice of program. Mail it to the 4050 Ser- 
ies Applications Library serving you; Library 
addresses are listed at the back of each 
TEKniques issue. 



51/00-0501/0 

51/00-0901/0 

51/00-1403/0 

51/00-1603/0 

51/00-4002/0 

51/00-5204/0 



51/00-6002/0 

51/00-8004/0 

51/00-8017/0 

51/00-8022/0 

51/00-9507/0 

51/00-9533/0 



* Documentation and listing only. 
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4050 Series Applications Library Program Abstracts 



Order 

Documentation and program listings of each program are 
available for a nominal charge. Programs will be put on tape or 
disc for a small recording fee per program plus the charge for the 
tape cartridge or flexible disc. One tape /disc will hold several 
programs. Programs will be recorded on like media only, i.e., 
programs on tape cannot be sent on disc and vice versa unless so 
noted in the abstract. 

(The program material contained herein is supplied without 
warranty or representation of any kind. Tektronix, Inc. assumes 
no responsibility and shall have no liability, consequential or 
otherwise, of any kind arising from the use of this program 
material or any part thereof. ) 

Domestic U.S. Prices: 

Documentation and listings $20 per program 

Recording Fee 5 per program 

Tape Cartridge 30 per tape 

Flexible Disc 15 per disc 



Contribute 

Contribute one program to the Library and receive three in 
exchange. Send in the membership card from your 4050 Series 
Graphic System Reference Manual to get the details. Or call us 
(503) 682-3411, ext. 3618. 

Forms 

Please use the Applications Library Order Form. Order forms 
are included in the Membership Packet and are available from 
your local Tektronix Sales Engineer. 

Outside U.S. 

Program contributions or orders outside the U.S. must be 
processed through the local Tektronix sales office or sent to one 
of the Libraries serving your area. See Library Addresses section. 



ABSTRACT #: 51/00-1407/0 

Title: Thevenin’s Theorem 
Author: John G. Lingle 
Tektronix, Inc. 

Beaverton, OR 
Memory Requirement: 16K 
Statements: 378 
Files: 1 ASCII Program 

The program provides computer aided 
instruction on how to THEVENIZE. 
The user inputs the number of voltage 
sources; the program randomly selects 
values for the voltages and resistors. The 
following may be selected to solve the 
problem: 

1. The program will progress through 
the solution, detailing each step. 

2. The user may solve the problem. 
The user has the option of using the 
4050 as a calculator to assist in com- 
puting the answer. Two attempts are 
allowed to input the correct answer. 
If the correct answer is not input the 
user may select to view the solution 
using the above process. 



ASSUME THE CIRCUIT IS OPEN WHERE 
THE "X" IS AND ONLY CONSIDER THE 
PORTION TO THE LEFT OF THE "X". 

TREAT THE CIRCUIT AS A SERIES BRANCH 
AND CALCULATE THE U0LTAGE AT THE 
POINT WHERE THE CIRCUIT WAS OPENED. 

PUSH RETURN: 



PUSH RETURN: 



01=330 03*220 

^Rl*203 ^R3*125 

-X-* — 



^R2«10 

02=-355 



THE EQUIMALENT RESISTOR IS 9.5 AND IS CALLED Rth (Resistance Thevenin). 

TO COMPLETE THIS STEP REPLACE THE TWO RESISTOR BRANCHES 
WITH THE PARALLEL EQUI0ALENT AND USE THE CALCULATED Ooc 
FOR THE 00LTAGE SUPPLY. 



PUSH RETURN: 



THERE IS JUST ONE RESISTOR BRANCH 
AND POWER SUPPLY LEFT THAT HAS NOT 
BEEN INCLUDED YET. 



PUSH RETURN: 



THIS IS THE FINAL Rth AND IS 8.8 



WOULD YOU LIKE TO SEE A THE0ENIH PROBLEM SOLOED <Y OR N>?Y 



ENTER NUMBER OF SUPPLIES 3 



THEOENIN'S is a nethod of solving circuits with 
TWO or nore power supplies. 



Using THEOENIN'S requires that only TWO supplies 
be considered at a tine. 



Each coFiplete step reduces the nunber of supplies 
and resistive branches by one. 



Eventually the circuit will be reduced to one 
equivalent supply and resistor. 

PUSH RETURN! 
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ABSTRACT #: 51/00-6122/0 3. Software communication protocol. 

Title: Data Communications Utility 4. Transmission error detection. 

Author: Chris Rotvik 

Tektronix, Inc. 

Irvine, CA 

Memory Requirement: 8K 
Peripherals: Opt 1 (RS-232) 

Statements: 135 
Files: 1 ASCII Program 

The program provides a means by which 
to transfer data files (or programs) be- 
tween 4050 Series Desktop Computers 
utilizing telephone lines and the optional 
Data Communications Interface (Opt 1). 

Features Include: 

1. Input and output files may be ASCII 
or binary. 

2. Interactive program allows operator 
selection of: 

A. Baud rate 

B . Transmit or receive 

C. Input or output file number 

D. Output file format (ASCII or bi- 
nary) if file format is not cur- 
rently defined 



ABSTRACT #: 51/00-6123/0 

Title: CYBER 175/UNIVAC 1108 Termi- 
nal Control 

Author: Chris Rotvik 

Tektronix, Inc. 

Irvine, CA 

Memory Requirement: 8K 
Peripherals: Opt 1 (RS-232) 

R06 Editor ROM Pack 
Statements: 63 each 
Files: 1 ASCII Program (CYBER 175) 

1 ASCII Program (UNIVAC 1108) 

The program consists of two separate 
programs. The first is CYBER 175 Ter- 
minal Control, the second is UNIVAC 
1108 Terminal Control. Each sets the 
proper optional Data Communications 
Interface (Opt 1) parameters and modes 
to interface a 4050 Series Desktop Com- 
puter to the specified host computer. 

Features are: 



C. Data send mode (to the system 
editor) 

D. Editor mode (utilizing the R06 
Editor ROM Pack) 

3. Designed around the need to per- 
form editing offline from the host 
computer. 



1 . Function key driven 

2. Provides the following modes of 
operation: 

A. Terminal mode 

B. Data receive mode (from the 
system editor) 
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ABSTRACT #: 51/07-8042/0 

Title: Coin Collection Inventory 
Author: Peter R. Hulick 
Lynchburg, VA 
Memory Requirement: 32K 
Peripherals: 4907 File Manager 
Statements: 354 
Files: 1 ASCII Program 

The program is used to create and main- 
tain an inventory for a coin collection. A 
file is created for each series of coins, 
one record contains each coin. The first 
record of each series contains the name 
of the series, denomination and number 
of items. Files may be reviewed, updat- 
ed, corrected and tabulated. 

The following options are available: 

1. Review of a series of coins 

Calls up a file of coins, displays the 
contents, tabulates and prints out 
the coins’ values. 

2. Add a new series of coins 

35 series of coins and suggested file 
identifiers are assigned in the pro- 
gram. These may be modified, add- 
ed to or deleted. 

3 . Change an existing series of coins 



PROGRAM TO INVENTORY COIN COLLECTION 

SELECT ONE OF THE FOLLOWING AND HIT RETURN KEY 

1 REVIEW OF A SERIES OF COINS 

2 ADD A NEW SERIES OF COINS 

3 CHANGE AN EXISTING SERIES OF COINS 

4 UPDATE THE PRICES OF AN EXISTING SERIES 

5 SUM THE VALUES OF ALL SERIES CURRENTLY ON DISK 

6 ADD TO A CURRENT SERIES 



YOUR CHOICE: _ 



»•* 


SERIES NAME 


FILE IDENTIFIER 


1 


TEST SERIES OF DIMES 


TEST10 


2 


MISC HALF CENTS 


MISHAFCENT 


3 


MISC LARGE CENTS 


MISLRGCENT 


4 


MISC 3 CENT PIECES 


M I SC3CENT 


5 


MISC 5 CENT PIECES 


MISC5CENT 


6 


MISC HALF DIMES 


MISHAFDIME 


7 


MISC DIMES 


MISDIMES 


8 


MISC 28 CENT PIECES 


MIS20CENT 


9 


MISC 25 CENT PIECES 


MIS25CENT 


10 


MISC 50 CENT PIECES 


MIS50CENT 


11 


MISC SILVER DOLLARS 


MISSILVD0L 


12 


US $1 GOLD 


USG0LD1 


13 


US $2-1/2 GOLD 


USG0LD2HAF 


14 


US $3 GOLD 


USG0LD3 


15 


US $5 GOLD 


USG0LD5 


16 


US $10 GOLD 


USGOLD10 


17 


US $20 GOLD 


USGOLD20 


18 


MEXICO 50 PESO GOLD 


MEXGOLD50P 


19 


MEXICO 2 PESO GOLD 


MEXG0LD2P 


20 


CANADA SOVEREIGN GOLD 


CANG0LDS0V 


21 


MISC SMALL CENTS 


MISCENTS 


22 


BUFFALO NICKELS 


BUFF5 


23 


JEFFERSON NICKELS 


JEFF5 


24 


BARBER DIMES 


BARBER10 


25 


MERCURY DIMES 


MERCURY10 


26 


ROOSEVELT DIMES 


ROOS10 


27 


BARBER QUARTERS 


BARBER25 


28 


STANDING LIB. QUARTERS 


STGLIB25 


29 


WASHINGTON QUARTERS 


WASH25 


36 


BARBER HALVES 


BARBER59 


31 


WALKIHG LIBERTY HALVES 


WALKING50 


32 


FRANKLIN HALVES 


FRANKLIN50 


33 


KENHEDY HALVES 


KENNEDY58 


34 


MORGAN SILVER DOLLARS 


MORGAN100 


35 


PEACE DOLLARS 


PEACE100 


ENTER FILE IDENTIFIER DESIRED 





Allows the correction of a coin in a 
series. 



4. Update the prices of an existing ser- 
ies 

Allows periodic updates of values of 
coins in a series. 



5. Sum the values of all series cur- 
rently on disk 

Lists all current series, each series 
total value, and the sum total of all 
series’ values. 

6. Add to a current series 



Allows the user to lengthen a series 
of coins that has been previously 
created. 



About 3000 coins may be stored on a 
disk. 



Tekniques 
Vol. 4 , No. 5 



21 






ABSTRACT #: 51/00-9534/0 

Title: Drafting Digitizer 
Author: Tom Sutherlin 

Cameron University 
Lawton, OK 

Memory Requirement: 16K 
Peripherals: 4662 Plotter 

Optional — 4956 Tablet 
Statements: 406 
Files: 2 ASCII Program 

The program consists of two parts; the 
menu and the digitizer. 

The menu is designed to be plotted on 
the 4662 Plotter. This is then placed on 
the 4956 Tablet to be used by the digitiz- 
er. The menu allows a blank area for 
sketching and digitizing with the basic 
symbols shown. 

The digitizer allows the user to transform 
sketches into a finished drawing using 
the basic symbols selected from the 
menu. Input is from the 4956 Tablet, but 
may be easily modified to allow input 
from the 4662 Plotter. The User-Defin- 
able Keys are used to select the screen 
or the plotter for output. 

200 X and Y coordinates may be input. 
For machines with more than 16K mem- 
ory the program may be modified to ac- 
cept additional X,Y coordinates. 

Menu selections include: 

Line — input points A and B, a line will 
be drawn between the two points 

Hidden Line — same as Line, but draws 
a dashed line 

Center Line — input points A and B, a 
dash .15 inches long will be drawn at the 
center distance of the line 

Circle — input point A (center) and B 
(right of center), used as radius to com- 
pute and draw circle 



Hidden Circle — same as Circle, but 
draws dashed circle 

Partial Arc — input points A, B, and C, 
calculates distance between A and B for 
radius, length is angle calculated be- 
tween AB and BC 

Partial Hidden Arc — same as Partial 
Arc, but draws a dashed line arc 

Leader — input points A and B, draws 
arrowhead at point A .07 inches wide 
and . 1 inches long with line to B 

Dimension Input — input points A and B 
(to the right), enter dimension from the 
keyboard (up to 16 alphanumeric charac- 
ters) data is output in a horizontal posi- 
tion 



End of Program — flashed “PROGRAM 
TERMINATED” on the screen and ends 
digitizing input 

Rectangle — input points A (lower left 
corner) and B (upper right corner), cal- 
culates the length of the horizontal and 
vertical sides 

Point-To-Point — input points A, B, C, 
D . . . , will draw straight line segments 
to connect the points, points may be at 
any position JmJ) 
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4050 Series Applications Libraries 



Africa, Europe, Middle East 

Contact local sales office 

Australia 

4050 Series Applications Library 
Tektronix Australia Pty. Limited 
Sydney 

80 Waterloo Road 
North Ryde, N.S.W. 2113 



Canada 

4050 Series Applications Library 
Tektronix Canada Ltd. 

P.O. Box 6500 
Barrie, Ontario 
Canada L4M 4V3 

Caribbean, Latin America 
and Far East (excl. Japan) 

IDD Group 
Export Marketing 
Tektronix, Inc. 

P.O. Box 500 
Beaverton, Oregon 97077 
U.S.A. 



Japan 

4050 Series Applications Library 
Sony/ Tektronix Corporation 
9-31 Kitashinagawa-5 
Tokyo 141 Japan 

United States 

4050 Series Applications Library 
Tektronix, Inc. 

Group 451 
P.O. Box 500 
Beaverton, Oregon 97077 
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